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ABSTRACT 

We report the discovery of a galaxy cluster at z=1.62 located in the Spitzer Wide-Area Infrared Extragalactic 
survey XMM-LSS field. This structure was selected solely as an overdensity of galaxies with red SpitzerllKAC 
colors, satisfying ([3.6]- [4.5])ab > ~0.1 mag. Photometric redshifts derived from Subaru XMM Deep Sur- 
vey (BV/z-bands), UKIRT Infrared Deep Survey-Ultra-Deep Survey (UKIDSS-UDS, y/T-bands), and from the 
Spitzer Public UDS survey (3.6-8.0 /im) show that this cluster corresponds to a surface density of galaxies at 
z « L6 that is > 20cr above the mean at this redshift. We obtained optical spectroscopic observations of galax- 
ies in the cluster region using IMACS on the Magellan telescope. We measured redshifts for seven galaxies 
in the range z=L62-L63 within 2.8 arcmin (< L4 Mpc) of the astrometric center of the cluster A posteriori 
analysis of the XMM data in this field reveal a weak (4(7) detection in the [0.5-2 keV] band compatible with 
the expected thermal emission from such a cluster The color-magnitude diagram of the galaxies in this clus- 
ter shows a prominent red-sequence, dominated by a population of red galaxies with (z-J) > 1.7 mag. The 
photometric redshift probability distributions for the red galaxies are strongly peaked at z = 1.62, coincident 
with the spectroscopically confirmed galaxies. The rest-frame (U-B) color and scatter of galaxies on the red- 
sequence are consistent with a mean luminosity-weighted age of L2 ±0.1 Gyr, yielding a formation redshift 
Z/ = 2.35 ± 0. 10, and corresponding to the last significant star-formation period in these galaxies. 

Subject headings: large-scale structure of the universe — galaxies: clusters: general — galaxies: clusters: 
individual (CIG J0218.3-0510) — galaxies: evolution 



1. INTRODUCTION 

Galaxy clusters provide important samples to study both the 
evolution of large-scale structure and the formation of galax- 
ies. The evolution of the number density of massive galaxy 
clusters involves primarily gravitational physics, which de- 
pends strongly on the cosmic mass density, the normalization 



' This work is based in part on observations made with the Spitzer 
Space Telescope, which is operated by the Jet Propulsion laboratory, Cal- 
ifornia Institute of Technology, under NASA contract 1407. This paper 
also includes data gathered with the 6.5 meter Magellan Telescopes lo- 
cated at Las Campanas Observatory, Chile. This work is based in part on 
data collected at Subaru Telescope, which is operated by the National As- 
honomical Observatory of Japan. 

^ George P. and Cynthia Woods Mitchell Institute for Fundamen- 
tal Physics and Astronomy, and Department of Physics and Astron- 
omy, Texas A&M University, College Station, TX, 77843-4242; pa- 
povich@physics.tamu.edu 

' Observatories, Carnegie Institution of Washington, 813 Santa Barbara 
St., Pasadena, CA, 91101 

Steward Observatory, University of Arizona, 933 N. CheiTy Ave., Tuc- 
son, AZ 85721 

^ W. M. Keck Postdoctoral Fellow, Harvard-Smithsonian Center for As- 
n-ophysics, 60 Garden St., Cambridge, MA 02138 

* Institute for Astronomy, Royal Observatory, University of Edinburgh, 
UK 

' Astronomy Centre, University of Sussex, Falmer. Brighton, UK 
Pontificia Universidad Catolica, Departamento de Astronomia y As- 
trofi'sica, 4860 Av. Vicufia Mackenna, Casilla 306, Santiago 22 Chile 

' Shanghai Key Lab for Astrophysics, Shanghai Normal University, 
Shanghai 200234, China 

Leo Goldberg Fellow, National Optical Astronomy Observatories, 
950 N. Cheny Ave., Tucson, AZ 85719 

" Department of Physics and Astronomy, University of Kansas, 1251 
Wescoe Hall Dr., Lawrence, KS, 66045-7582 

Argelander Institute for Astronomy, Bonn University, Auf dem HUgel 
71, 53121 Bonn, Germany 

" Service d'Astrophysique, CEA Saclay, 91191 Gif sur Yvette, France 



and shape of the initial power s pectrum, as well as on the dark 
energy equation of state (e.g.. lEke et al. [19981 iBa hcall et al 



1999; Borgani et al.ll2001l:lHaiman et all 2001l:ISpringel et 
20051; iPacaud et al.lf2007riJee et alj|2009l) . Because the mas 
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sive galaxies in clusters formed nearly contemporaneously 
at z :: j> 1 with similar star-formation and a s sembly histories 
(e^ lStanford et al J [19981 iTran et al J [20071 lEisenhardt etll] 
120081) . observations of the evolution of distant cluster galax- 
ies provide strong constraints on hierarchical galaxy evolution 
models, which ma ke detailed predictions for the formation of 
these objects (e.g.. lDe Lucia & Blaizotl[2b07h . 

Studies of high-redshift clusters have been frustrated by 
small sample sizes. Few clusters have confirmed redshifts 
beyond z - 1.3 (e.g. iMulUs et al.|[2q()5l; [Stanford et al. ll2005l; 
Brodwin et al] 120061 [Stanford et alj 1200 6; Eisenhar dt et alj 
20081 iKurk et al.1 120091: [Wilson et al.1 120091) . This dearth 
of detected galaxy clusters at z > 1.3 stems from the sig- 
nificant challenges and potential biases in identifying these 
structures. Deep X-ray surveys have identi fied spectro- 
scopi c ally confirmed cluste rs to z ^ 1.5 (e.g., iRosati et"an 
[20041; [Stanford et"al] 120061) . but X-ray selection typically 
require relaxed systems, which is unlikely to be the case 
at high redshifts where cluster progenitors will be less 
massive and more disordered, with less time for the in- 
traclu ster medium (ICM) to thermalize (e.g., iRosati et"an 
120021) . Searches for galaxy overdensities around distant ra- 
dio g alaxies requ ire the presence of a mass ive central galaxy 



je.g., Kurketal. 2000; Milev et al. 2004; Stern et al. 2003 
| Kodama et al. 2007: Venemans et al. 2007; Zirm et al, 2008; 
Chiaberge et al which may not be an intrinsic prop- 



erty of cluster progenitors 

Other searches utilize the 
color-magnitude relation 



empirically observed, tight 
in central cluster galaxies 
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(the "red sequence" e.g.. IVisvanathan & Sandag^ 119771: 
Gladders & Ye^ llOOl iGladders et al.l 120071: iKaiisawa'etal] 
20061: iMuzzin et al.l 120091) . However, this selection is 



biased potentially against clusters whose galaxies have 
had recent star formation (and overdensities of blue, star- 
form ing galaxies at h igh redshift have been identified, 
e.g., ISteidel et al.l l2005h . Most studies of cluster galaxies 
imply their stellar populations form ed at z/ > 1.5 (e.g., 
Ivan Dokkum & van der M^ell 120071) . concurrent with the 
peak epoch in the star-formation rate density from UV and 
IR measurements (see e.g., Hopkins & Beacom 2006). As 
searches for clusters approach the redshift of their formation, 
cluster galaxies should show increasing indications for star 
formation with less time available to build up a substantial 
red-sequence population. 

We have initiated a search for galaxy cluster candidates at 
z > 1.3 selected solely as overdensities of galaxies with red 



colors using dat a from the Infrared Arra y Cam- 
Fazio et aLll2004l) o n board Sp itzer ( Wern er et al.l 



[3.6] -[4.5 

era (I RAC, 

120041) following the method of lPapovichI (12008.) PI At z<l 
model stellar populations have blue [3.6] -[4.5] colors be- 
cause these bands probe the stellar Rayleigh-Jeans tail (with 
the expection of some IR-luminous, star-forming galaxies at 
z ~ 0.3 which have a contribu tion of warm dust to their near- 
IR colors, see lPapovichl2008h . At z> 1 both star-forming and 
passively evolving stellar populations appear red in [3.6]- 
[4.5] as these bands probe the peak of the s tella r emission 
at 1.6 nm (see lSimpson&Eisenhard3ll999l: ISawicki 2002; 
iPapoviclij 12008*). Therefore, selecting overdensities of red 
[3.6]- [4.5] sources potentially identifies high-redshift clus- 
ter candidates wi th little bias from galaxy stellar populations. 
iPapovichI (12008!) showed that IRAC-selected z > 1.3 clus- 
ter candidates from the Spitzer Wide-Infrared Extragalactic 
(SWIRE) survey have clustering scale lengths of ro ~ 20 /i"' 
Mpc, consistent with other high-redshift galaxy clusters (see 
iBrodwin et al]|2007b . 

Here we report the discovery and spectroscopic confirma- 
tion of a galaxy cluster, CIG J0218. 3-0510, (corresponding to 
Infrared Cluster "A", 1R C-0218A, in the 0218-051 field, se- 
lected in lPa povichll2008l) , and we discuss its photometric and 
spectroscopic properties. This galaxy cluster was identified 
using our IRAC color selection with no other additional crite- 
ria imposed. In § 2 we discuss the target selection, and pho- 
tometric and spectroscopic observations. In § 3 we present 
the evidence supporting the assertion that this structure is a 
galaxy cluster, and in § 4 we discuss its properties. In § 5 
we summarize our results. Unles s otherwise note d we report 
all magnitudes in reference to the iJohnsonI (119661) magnitude 
system relative to Vega. We explicitly de note magnitudes rel- 
ative to the absolute bolometric system (lOke & Gunnlll983h 
with an AB subscript, m^B = 23.9-2.51og(/,y/l /iJy). We use 
cosmological parameters 17,,, = 0.3, A = 0.7, and H = 10 km 
s~' Mpc~' throughout. For this cosmology, the angular diam- 
eter distance is 0.5 Mpc arcmin"' at z = 1 .62. 

2. SAMPLE SELECTION AND SPECTROSCOPIC OBSERVATIONS 

2.1. Identification of High-redshift Galaxy Cluster 
Candidates 

We identified galaxy cluster candidates at z > 1.3 using 
the IRAC data from the SWIRE survey. These data cover 

Throughout we denote magnitudes measured in the 3.6, 4.5, 5.8, and 
7.9 fim IRAC channels as [3.6], [4.5], [5.8], and [8.0], respectively. 



roughly 50 deg^ divid ed over six fields separated on the sky 
dLonsdale et al.ll2003l) . We used a simple color selection to 
identify high-redshift galaxies fro m Spitzer/IRAC 3.6 and 
4.5 /im photometry. As discussed in iPapovicE (120081) >90% 
of galaxies with z > 1.3 have ([3.6] - [4.5])ab > -0.1 mag. 
[Papovich (2008) identified candidate galaxy clusters at z > 1.3 
by selecting overdensities of > 30 objects satisfying this color 
criterion within radii of I'A (corresponding to r < 0.7 Mpc 
at z = 1.5) from the SWIRE surv ey data. CIG J0218.3-0510 
was identified in IPapovichI (l2008h in the SWIRE XMM-LSS 
field, and has astrometric coordinates, a = 2''18'"21.3', S = 
-05°10'27" (J2000), derived from the centroid of the SWIRE 
IRAC sources in this overdensity. The left panel of figure [T] 
shows a false-color image using the B-band (blue), /'-band 
(green) and 4.5 fim image (red) of the field centered on the 
coordinates of the IRAC-selected overdensity. 

2.2. Optical and Near-IR Imaging 

The XMM-LSS includes deep optical (0.4-1 pm) and near- 
IR imaging (1-2 pm) in a portion of the field covering 
roughly 0.70 deg^. Optical imaging is available from th e 
Subaru-XMM Deep Survey (SXDS, iFurusawa et"aD 120081) . 
obtained with the Subaru Prime Focus Camera (Suprime- 
Cam) in the broad bandpasses B, R, ;', and z'. Near-IR imag- 
ing is available from the UKIRT Infrared Deep Sky Survey 
(UKIDSS, data release 1, Lawrence et al. 2007) ultra deep 
survey (UDS) in the broad bandpasses J and K. For the 
work here, we utiliz ed the jiT-selected cata log of the SXDS 
and UDS data from IWilliams et al.l (l2009h . These catalogs 
reach 5cr-limiting magnitudes in 1 "75-diameter apertures of 
Bab < 27.7,/? AB < 27.1, /AB < 26.8, zab < 25.5, 7ab < 23.9, 
and Kab < 23.6 mag. The K-selected catalogs also include 
quasi-total K magnitudes, measured in elliptical apertures 
based on the hght profile for each so urce. For details, we refer 
the reader to IWilliams et al.l (l2009h . We corrected the fixed- 
aperture magnitudes to quasi-total magnitudes using a unique 
aperture correction for each source, defined as the difference 
between the fixed aperture and total magnitude derived from 
the /T-band for that source, K{ap)-K{tot). We then applied 
the aperture correction for each source to its fixed aperture 
magnitudes measured in the other SXDF and UKIDSS bands. 

2.3. Spitzer IRAC Imaging 

While the SWIRE IRAC data are sufficient for the selec- 
tion and study of the galaxies in this cluster, we made use 
of deeper IRAC imaging in this field from the Spitzer public 
legacy survey of the UKIDSS UDS (SpUDS, PI: J. Dunlop). 
The SpUDS data cover a field of area 1.5 deg^ to deeper lim- 
iting flux densities than those available with the SWIRE sur- 
vey data. We measured photometry in 4"-diameter apertures 
in each of the IRAC 3.6, 4.5, 5.8, and 8.0 ^m images. We 
applied correction factors of 0.32, 0.36, 0.55, 0.68 mag, re- 
spectively, which corrects these aperture magnitudes to total 
magnitudes for point sources. We matched sources in each 
of the IRAC-selected catalogs to those in the UKIDSS UDS 
/T-selected catalog within 1" radii. We combined the "to- 
tal" aperture-corrected magnitudes for the IRAC data with the 
aperture corrected photometry for the SXDS and UDS data. 

2.4. Photometric Redshifts and the Integrated Redshift 
Probability 

The merged SXDS, UDS, and IRAC catalog covers a wave- 
length baseline of 0.4-8 fj,m, and we used these data to 
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Figure 1. The left panel shows a false-color image of the tai'get field. Blue corresponds to the Suprime-Cam B-band, Green to the Suprime-Cam j-band. 
and Red to the Spitzer 4.5 band. The images have not been corrected for variations in the data image quality. The image spans 3' X 3', corresponding to 
1.5Mpc X 1.5Mpc at 2 = 1.62. The heavy white bai' shows a distance of 0.7 Mpc at the redshift of the cluster. Small green circles denote candidate cluster 
members with > 0.5 as defined in § 12. 4| Orange squares denote those objects with spectroscopically confirmed redshifts 1.62 < z< 1.65. The contours denote 
regions with 5, 10, and 15(t above the mean density of galaxies with 1.5 < Zphot < l -^- The right panel shows the surface density of galaxies with 1.5 < jphot < 1-7 
in units of the number of standard deviations (cr) above the mean density, ranging from -1 to 20 (as indicated in the plot legend). 



Study the cluster-candidates selected from the SWIRE IRAC 
data. We derived photometric redshi fts for each source i n 
the A'-selected catalog using EAZY (iBrammeret al.llIOOSh . 
We used the default galaxy spectral energy distribution tem- 
plates with a /T-band prior based on the luminosity func- 
tions of galaxies in a semi-analytic simulation. We de- 
rived the most likely photometric redshift as well as the full 
photometric-redshift probability distribution function, Piz), 
normalized such that J P{z)dz = 1 when integrated over all 
redshifts. Our comp arisons against the spectroscopic redshifts 
in the SpUDS field (lYamada et aU lMBEr Simpson et alM M. 



Ivan Breukelen et alj|2007ll2009l C. Simpson, in preparation) 
showed that the most likely photometric redshifts have un- 
certainties of A([zsp-Zph]/[1 +Zsp]) = 0.04 derived froin the 
normalized median absolute deviation (iBeers et al.lll990h for 
the more than 200 galaxies (excluding broad-line AGN) with 
spectroscopic redshifts in the range 1 .0 < z < 2.0. 

To increase the efficiency of our spectroscopic observa- 
tions, we computed the surface density of galaxies in the 
SpUDS field in coarse redshift intervals, and we prioritized 
those IRAC-selected cluster candidates (selected over the 
much larger SWIRE field) that corresponded also to large 
overdensities in photometric redshift. To measure the surface 
density of galaxies, we divided galaxies into redshift intervals, 
Az = 0.2, and measured the angular distance from each object 
to the seventh-nearest neighbor, dj, and then computed the 
corresponding surface density, cx {d-j)~^. We tested other 
definitions for the nearest neighbor, which produced similar 
results (changing the definition of the Mh-nearest neighbor 
changed primarily the angular resolution of the surface den- 
sity map). We then calculated the mean and standard devia- 
tion of the surface density across the entire UDS field. 

CIG J0218. 3-0510 appears as a strong overdensity of galax- 
ies with 1.5 < Zph < 1-7. Figure [T] (right panel) shows the 



surface density of galaxies in this photometric-redshift range 
centered on the IRAC-selected overdensity CIG J0218.3- 
0510. The color shading corresponds to the number of stan- 
dard deviations above the mean surface density of galaxies 
at this redshift over the UDS field. CIG J0218.3-0510 corre- 
sponds to a > 20(7 surface density of galaxies in this redshift 
interval. 

Many of the galaxies in the field of CIG J0218.3-0510have 
photometric redshifts centered tightly around Zph ~ 1.6. We 
quantified the likelihood of galaxies being associated in red- 
shift by defining the integrated redshift probability. 



V,= / P(z)dz, 



(1) 



integrated over z = Zcen ± &■ For CIG J0218.3-0510, we 
used Zcen = 1-625 and 5z = 0.05 x (1 +Zcen) (therefore in- 
tegrating over 1.49 < z < 1.76), approximately the 68% 
confidence range on the photometric redshifts. This was 
motivated by other cluster-member selec tion methods us 



ing photometric-redshift-selected samples ( Brunner & Lubir 



20001: iHalUdav et al.ll2004lEisenhardt et al 



l2008tlPell6et"ar 



20091) . Simply defined, the integrated redshift probabiUty is 



the fraction of the photometric redshift probability distribu- 
tion function within th e redshift intervals (see discussion in 
iFinkelstein et alJl2010h . An integrated probability of = 0.5 
means that 50% of the integrated photometric redshift distri- 
bution lies between Zcen -5z < z < Zcen + 5z. The galaxies in 
the field of CIG J0218.3-0510 with V:, > 0.5 are indicated in 
figure [U 

2.5. Spectroscopic Observations 

We targeted galaxies with high integrated redshift proba- 
biUty in the region of CIG J0218. 3-0510 using the Inamori 
Magellan Areal Camera and Spectrograph (IMACS) on the 
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Table 1 

Spectroscopic Redshifts in tlie CIG J0218.3-0510 Field 



R.A. 


Decl. 




z 


cr- 


R 


A 


r 


(J2000) 


(J2000) 








(mag) 


(') 


(Mpc) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


2''18'"22.30' 


-5°10'34.7" 


0.50 


1.6224 


0.0005 


22.9 


0.28 


0.14 


2''18'"20.37' 


-5°09'56.3" 


0.42 


1.6303 


0.0012 


23.6 


0.57 


0.29 


2''18'"19.16' 


-5° 10' 14.7" 


0.47 


1.6230 


0.0005 


23.0 


0.57 


0.29 


2''18"'24.14' 


-5°09'45.2" 


0.46 


1.5356 


0.0004 


23.2 


0.99 


0.51 


2''18'"24.48' 


-5°09'39.7" 


0.40 


1.6228 


0.0005 


23.4 


1.12 


0.57 


2* 18™ 17.20' 


-5°11'05.7" 


0.61 


1.6487 


0.0006 


23.8 


1.21 


0.61 


2''18'"23.92' 


-5°09'20.6" 


0.54 


1.6222 


0.0005 


22.8 


1.29 


0.65 


2'' 18'" 18.92' 


-5°08'00.3" 


0.32 


1.6234 


0.0013 


23.7 


2.52 


1.28 


2'' 18'" 17.46' 


-5°08'02.2" 


0.48 


1.4962 


0.0003 


22.9 


2.60 


1.32 


2* 18'" 15. 18' 


-5°08'11.9" 


0.64 


1.6224 


0.0011 


23.8 


2.72 


1.38 


2'' 18'" 14.50' 


-5° 06' 59.1" 


0.47 


1.6094 


0.0008 


23.5 


3.86 


1.96 



Note. — (1) Right ascension, (2) Declination, (3) Integrated redshift probability, 
(4) spectroscopic redshift, (5) redshift uncertainty, (6) Suprime-Cam S-band mag- 
nitude, (7) angular separation between the galaxy and the cluster astrometric center 
(see § 2.1), (8) projected physical separation between galaxy and the cluster center 
for z= 1.62. 



Magellan/Baade 6.5 m telescope on 2008 Oct 30-31, 2008 
Nov 18-19, and again on 2009 Sep 11-12. At the f/2 focus 
IMACS provides multiobject spectroscopic observations over 
a 27.2'-diameter field of view (FOV), which allowed us to tar- 
get roughly 100 galaxies in as many as 10 cluster candidates 
using a single slitmask. We targeted two separate IMACS 
fields within the UDS, covering 10 cluster candidates in field 
1 and 7 in field 2, using two slitmasks per field to alleviate slit 
collisions. CIG J0218. 3-0510 was one of our highest priority 
targets given the high surface density of objects with high V^- 
The other cluster candidates will be discussed in a forthcom- 
ing paper. 

We observed with IMACS using the 200 lines/mm grat- 
ing with the OG570 blocking filter, which provided ^7 A 
resolution covering 0.6 - 1 /im. The IMACS CCDs were 
upgraded in 2008 and have very high red sensitivity, with 
about a factor of two improvement in sensitivity beyond 
8500 A (A. Dressier 2008, p rivate communication). We used 
the "n od-and-shuffle" mode ( iGlazebrook & Bland-HawthornI 
12001b with 1" X 2.2" slitlets, which greatly improves the 
background subtraction and facilitated our spectroscopic red- 
shift success, especially for galaxies at z > 1.3. For practical 
purposes, we prioritized galaxies with 7? < 23.3 mag, but in- 
cluded galaxies toR< 23.8 mag (and a few even fainter galax- 
ies). Conditions each night were generally clear, although 
not photometric, and some time was lost to poor photometric 
conditions (including most of the time during the 2009 Sept 
run). Typical image quality during good conditions ranged 
over 0.6-1" at full width at half maximum for point sources. 
Individual exposure times were 1800 s. The total co-added 
exposure times of data taken under good photometric condi- 
tions varied for each mask, and were 3-4 hours on source. 

We reduced the IMACS spectroscopic data using the 
Carnegie Observatories System for MultiObject Spectroscopy 
(COSMOS, v2.13l3 The reduction steps to produce 2D 
spectra include wavelength calibration, bias subtraction, flat 
fielding, sky subtraction, co-adding the separate frames, and 
cosmic-ray removal. We extracted ID spectra at each nod- 
and-shuffle position separately, as well as the associated un- 
certainty on the spectra propagated through the reduction 

" http://obs.camegiescience.edu/Code/cosmos 



pipeline. We then coadded the individually extracted ID spec- 
tra for each target. We used spectrophotometric standards 
taken at the end of each night to provide flux calibration. Be- 
cause photometric conditions varied over the course of each 
night, our photometric calibration is not absolute (although 
this does not affect our ability to measure redshifts). 

We took spectra of 24 sources with 7? < 23.8 mag, inte- 
grated redshift probability Vz > 0.3, and with an angular sep- 
aration of A < 4' of the astrometric centroid of CIG J0218.3- 
0510 (corresponding to a physical separation of r < 2.0 Mpc 
at z = 1.62). We measured spectroscopic redshifts for 11 of 
these galaxies, a redshift success rate of about 45% (11/24). 
These redshifts are listed in table [T] In all cases, the redshifts 
are secured on the basis of the [O II] emission line. We in- 
spected the 2D reduced data and verified that the emission 
feature is present in both independent spectra obtained from 
the nod-and-shuffle observation. All of these galaxies were 
targeted as candidate members of this galaxy cluster Seven 
of the galaxies have 1 .62 < z < 1 .63 with a mean (z) = 1 .622. 
Of the remaining galaxies, two have z = 1 .649 and 1 .609 (ve- 
locity separations of <3100 km s~'), and the other two have 
z = 1.536 and z = 1.496. Including all nine galaxies within 
3100 km s"', we derive a mean redshift (z) = 1.625. 

Figure |2] shows the one-dimensional spectra for six galax- 
ies with 1.62 < z < 1.65 with the smallest angular separation 
(A < 1.'4) from the astrometric center of CIG J0218.3-0510. 
Each spectrum shows the presence of [O II] in emission. Sev- 
eral of the spectra show evidence for interstellar absorption 
from Mg II and possibly Fe II (one galaxy at z = 1 .649 shows 
weak Mg II emission, likely from an AGN). 

As the rest of the absorption features are weak, we show in 
figure [3] a weighted-mean, composite spectrum for the seven 
galaxies with 1 .62 < z < 1 .63, all within r < 1 .4 Mpc. To con- 
struct this spectrum, we scaled each spectrum to the mean flux 
density in the rest-frame wavelength range 2600-3100 A. We 
then weighted each spectrum by the inverse variance using the 
uncertainty spectrum. The stacked spectrum shows broad fea- 
tures from ^3200-3400 A, which are artefacts of the stack- 
ing procedure. Because we have scaled each spectrum to the 
mean flux, the fainter galaxies contribute more noise to the 
stack. These features correlate with strong sky emission, and 
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Figure 2. Extracted IMACS spectra for the six galaxies witli 1 .62 < z < 1 .65 witliin 1 .4' of the cluster center (coiTesponding to r < 0.7 Mpc at z = 1 .62). The 
inset panel shows the region around [O II] A3727 at the measured redshift. In all panels, the black line is the measured spectrum. Gaps in the spectra show 
wavelengths that fall on the "gaps" between the IMACS CCDs. The red line shows the Icr uncertainties propagated through the data reduction. The thick black 
lines show the expected locations of Fe II A2374, Mg II A28()0, and [O II] A3727 for the measured redshift. 

a result of the higher noise in this region of the spectrum. The 
stacked spectrum shows a strong upturn in flux density at the 
red end of the spectrum. This feature is unlikely to be intrin- 
sic to the galaxies, and arises from the fluxing uncertainties 
at the red end of the spectra, due to the decreasing sensitivity 
of the CCD and the red cut-off of the blocking filter Nev- 
ertheless, strong emission from [O II] A3727 is observed in 
the composite spectrum, which is a result of the spectroscopic 
identification process. The stacked spectrum shows several 
strong absorption features, in particular Mg II AA2976,2804, 
Mg I A2851 and absorption from various blends of Fe II lines. 

The narrow redshift range of the galaxies in table [T] sug- 
gests they are physically associated. Using the spectroscopic 
redshifts to infer the dynamical conditions of the cluster is 
dubious because of the unknown dynamical state (see § 4.2). 
If we assume the galaxies within <1500 km s"' of the mean 
redshift and within 0.9 Mpc of the astro me trie center of the 
IRAC-selected overdensity sample adequately a virialized 
structure, then their redshifts correspond to a line-of-sight 
veloci ty dispersion of a = 8 60 ± 490 km s"' using the defini- 
tion in lCarlberg et al.l (119961) . The uncertainty is derived from 
a bootstrap resampling of the data (and we make no attempt 
to remove the additional uncertainty from the redshift errors). 
However, we further qualify this velocity dispersion because 
of several biases that likely affect the redshift success rate 
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Figure 3. Composite one-dimensional IMACS spectra as a function of rest- 
frame wavelength for the seven galaxies within r< 1.4 Mpc of the astrometiic 
center of CIG J0218.3-0510 with spectroscopic redshifts 1.620 << z < 1.630 
(see Table [T). The black line shows the composite, weight-averaged spec- 
tram, and the red Une is the weighted uncertainty, dominated by emission 
from the sky. Strong [O II] is observed, as well as several absorption fea- 
tures, including the doublet Mg IlAA2976,2804, and absorption from blends 
of Fe II line. The absorption feature at rest-frame 2900 A corresponds to 
telluric absorption at »:i7600 A in the observed frame. 
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Figure 4. Summed XMM image centered on the position of CIG J0218.3- 
0510. CIG J0218.3-0510 is located on the edges of three adjacent XMM 
positions. The coadded XMM data yield an unambiguous source with 4cr 
significance in the [0.5-2.5] keV band within 10" of the position of the clus- 
ter. The circle indicates the associated X-ray source and has a diameter of 
40". The units of the color scale is erg s"' cm"^ arcmin"^ as indicated by the 
running bar along the bottom of the image. 

for galaxies in this cluster. Firstly, we prioritized only those 
galaxies with R < 23.8 mag for spectroscopy, preferentially 
excluding fainter galaxies. Second, even with this magnitude 
criterion, most of the galaxies are faint (23 <R <24 mag) and 
all the redshift identifications above are based primarily on the 
[O II] A3727 emission line. At z = 1 .62- 1 .63, this line falls at 
9780-9800 A, adjacent to a strong sky line at 9800 A, which 
hinders the identification of galaxies at z w 1 .63 with emis- 
sion lines (the only galaxy at this redshift identified above in- 
cludes strong Mg II absorption). Therefore, our spectroscopic 
redshifts will be biased against objects at z « 1.63. If the 
cluster lies at z = 1 .625 then roughly half the galaxies with 
redshifts in the upper end of the distribution will be prefer- 
entially missed, biasing the mean redshift and the velocity 
dispersion measurement. Therefore, while we conclude the 
galaxies are physically associated as a galaxy (proto-)cluster, 
we caution against too much interpretation of the dynamics of 
this cluster using the redshifts above. 

2.6. XMM-Newton X-ray Imaging 

CIG J0218.3-0510 is located within the 1° X MM Subaru 
Dee p Field enclosed in the XMM-LSS survey (iPierre et al.l 
12004). While our cluster selection method did not involve 
X-ray criteria, the presence of diffuse X-ray emission at the 
cluster location would provide independent confirmation of 
the existence of a deep g ravitational potent i al wel l. In a re- 
cent analysis of these data lFinoguenov et al.l (120101) identified 
two diffuse X-ray sources at distances of lf8 and 2. '9 from 
CIG J0218.3-0510, with redshift estimates of z = 1.6- 1.8. 
These may be associated with CIG J0218. 3-0510 although 
their large offset in angular distance challenges this interpre- 
tation. 

Following our discovery of the close association in redshift 
of the galaxies in CIG J0218. 3-0510, we performed a poste- 
riori analysis of the XMM data at the location of this object. 
CIG J0218. 3-0510 is present on three adjacent XMM point- 
ings, near the edge of each of them (off-axis values >12'). We 
coadded the three XMM images centered on the position of 
CIG J0218.3-0510, and applied an adaptive smoothing. The 
X-ray image is shown in figure]?] The image has an effective 



exposure time of w 3 x 12 ks at the position of CIG J0218.3- 
0510. There is an unambiguous source within a radius of 10" 
of the cluster position detected with a 4cr significance. Within 
a 40"-diameter aperture we measured 65 photons (±30%) in 
the [0.5-2 keV] band. Owing to the low number of photons 
and the large off-axis position of CIG J0218. 3-0510 in the 
XMM images, we are unable to constrain accurately the spa- 
tial extent of the X-ray source nor to estimate a gas tempera- 
ture. However, the detection of the X-ray emission provides 
an estimate for the virial mass of CIG J0218. 3-0510, which 
we discuss in § 4.2. 

3. THE NATURE OF CIG J0218.3-0510 

In this section we discuss the evidence that CIG J0218.3- 
0510 is a galaxy (proto-)cluster Strictly speaking, a "cluster" 
is an object that is fully virialized, while a "proto-cluster" is 
an object that will eventually virialize at later times (lower 
redshift). It remains to be seen if the dark-matter halo of 
CIG J02 18.3-05 10 is fully virialized or if it is still assembling. 
Nevertheless, we will blur the distinction between these defi- 
nitions and use the term "cluster" to mean both. 

The photometric redshift distributions and spectroscopic 
redshifts of galaxies in CIG J0218. 3-0510 indicate a large 
overdensity of galaxies at z = 1 .62 with a high surface density 
within r < 0.7 Mpc. This is illustrated in figurejl] The galax- 
ies with spectroscopic redshifts in this region of high surface 
density show a preponderance of sources around z = 1.62. 
Furthermore, the IRAC-selected galaxy members of this clus- 
ter have high integrated redshift probability as defined in 
equation 1, implying a high likelihood of being at the clus- 
ter redshift. Galaxies with Vr > 0.5 are indicated with circles 
in figure JT] This includes many red galaxies, which are con- 
centrated near the center of CIG J0218. 3-0510 and have very 
high Vz values. 

Figure]5]shows a montage of photometric redshift probabil- 
ity distribution functions, P(z), for the seven galaxies with the 
highest integrated redshift probabilities, V^, with an angular 
separation of A < 1' of the cluster center Without excep- 
tion these galaxies have red (z-J) colors. We targeted sev- 
eral of these galaxies with spectroscopy from IMACS during 
our 2009 observing run. However, their spectra show only 
faint optical continua, with no discernible emission features. 
Nevertheless, while we are unable to derive spectroscopic red- 
shifts for these galaxies, the lack of any emission features is 
consistent with their photometric redshifts. The photometric 
redshifts of these galaxies are driven by the strength of the 
apparent 4000 A/Balmer break redshifted between the z- and 
7-bands at z = 1 .62. As a result, the photometric redshift P{z) 
functions are sharply peaked at this redshift, and these galax- 
ies have Vz > 0.75, implying that more than 75% of their red- 
shift probability lies between 1.49 < z < 1.76 with a most 
likely redshift at z — 1 .62, the mean spectroscopic redshift. 

In contrast, the galaxies for which we obtained spectro- 
scopic redshifts have lower V^- These galaxies all show [O II] 
in their spectra, implying relatively recent star formation. 
Therefore they have weaker 4000 A/Balmer breaks and sub- 
sequently less-well constrained P(z). This is illustrated in fig- 
ure]5] which shows the P{z) for the seven galaxies with spec- 
troscopic redshifts at that of the cluster All of these galaxies 
have bluer (z-J) colors, consistent with the evidence of more 
recent star formation. 

The red galaxies in this overdensity dominate the color- 
magnitude relation for this cluster, simi lar to the relations 
observed in clusters at z < 1.4 (e.g., iBower etal] 119921 : 
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Figure 5. Photometric-redshift probability-distribution functions, P(i), for 
galaxies in CIG J0218. 3-0510. The left figure shows the P{z) for the seven 
galaxies that have integrated redshift probabilities, Vz, > 0.75 with angular' 
separations of A < 1 ' of the cluster center. Without exception, these galaxies 
all have very red (z-J) colors. Each panel gives the {z—J) color, 7-magnitude, 
the most likely photometric redshift, the integrated redshift probability, and 
angular separation from the cluster center. Although these galaxies lack spec- 
troscopic information, their large integrated redshift probabihties support the 
assertion that they are at the cluster redshift, z = 1.62, indicated by the ver- 
tical dashed lines. The cluster galaxies with spectroscopic redshifts all show 
[O II] in their spectra, implying ongoing star formation (see figure|2). They 
have bluer (z-J) colors and lower values of V-, presumably because they 
have weaker 4000 A/Balmer breaks owing to the star formation activity. The 
light figure shows the P(z) for the seven galaxies with spectroscopic redshifts 
1.62 < z < 1.63. Each panel shows the same information as for the left fig- 
ure, except that they give the spectroscopic redshift instead of the most likely 
photometric redshift. 



Ellis et al I [T99I [Stanford et all [T991: Ivan Dokkum et al I 
1998allbt iBlakeslee et al.l 12003'. '2006^, 'De Lucia et al.l 120071: 
Tran et al.l 120071: iLidman et al. 2008; Mei et al. 2009). Fig- 
ure |6] shows a z-J versus / color-magnitude diagram for all 
galaxies within 2 arcmin (r < 1 Mpc at z = 1 .62) of the astro- 
metric center of the IRAC-selected overdensity. The z and J 
bandpasses span the redshifted 4000 A/Balmer break at red- 
shift z= 1.62 providing strong contrast between the cluster 
galaxies and those in the field. Field galaxies are shown as 
small data points in figure|6] and these have "P- < 0.3, imply- 
ing they are foreground or background galaxies. Larger sym- 
bol sizes correspond to higher values of P^, implying those 
objects have a greater likelihood of being at the cluster red- 
shift. The cluster galaxies on the red-sequence typically have 
the highest values of P~, especially at the bright end where 
photometric uncertainties have a smaller effect on the colors. 
The well-defined color-magnitude relation in CIG J0218.3- 
0510 is evidence for the cluster-like nature of this object. 
This is the highest redshift spectroscopically confirmed clus- 
ter with such a well-defined red sequence. 



The red galaxies that dominate CIG J0218. 3-0510 are cen- 
trally concentrated, and nearly all are within A < 1' of the 
cluster center. Very few of the galaxies on the red sequence 
have angular separations >1'. Of the few that have larger 
angular separations, most reside within the surface-density 
contours that extend slightly to the north-east and north- 
west from the central region (see figure [T]). This property of 
CIG J0218.3-0510 is consistent with the color-density rela- 
tions observed i n other clusters at low and high redsh ifts (e.g., 
lDressleii[T980t iLidman et al']l2008t iMei et al.1 llool provid- 
ing further support for the cluster-like nature of CIG J0218.3- 
0510. 

In summary, the spectroscopic and photometric redshifts, 
the color-magnitude relations, and the surface density-color 
relations all support the conclusion that CIG J0218. 3-0510 is 
a galaxy cluster This interpretation is reinforced by the faint 
X-ray emission measured at the location of this object (see 
also § 4.2). We therefore conclude these galaxies are physi- 
cally associated with each other as a cluster, even though we 
do not have the data to determine if they are fully virialized. 
In the next section we discuss the properties of this object and 
the cluster galaxies. 

4. DISCUSSION 

From the evidence presented above, we conclude that 
CIG J0218. 3-0510 represents a galaxy cluster at z = 1.62. 
Although other candidates for high-redshift galaxy (proto- 
)clusters at z > 1.5 have been reported (e.g., [Milev et al. 
200?; Brodwin et al." 2007; McCarthv et al. 2007; Zirm et al 
2008; J Eisenhardt et a l. 2008^ Andreon et al. 2009; Kurk et al 
2009t IChiaberge eFall l2010l) . this is the highest redshift, 
spectroscopically-confirmed cluster with a strong, well- 
defined red sequence. While the definition of galaxy "cluster" 
may be subject to semantics (see above), CIG J0218. 3-0510 
shows all the characteristics indicative of lower-redshift, rich 
galaxy clusters. 

The fact that CIG J0218.3-0510 has a well-defined red se- 
quence of strongly clustered red galaxies is not a result of the 
selection method. CIG J0218. 3-0510 was selected as an over- 
density of galaxies with red IRAC colors (see § 2.1). The 
IRAC colors at this redshift are blind largely to variations of 
rest-frame optical colors, and should be sensitive to galax- 
ies with both re st-frame red and blue UV-optical colors (see 
lPapovichll2008l) . With our full sample, we will compare the 
properties of CIG J0218.3-0510 against other overdensities of 
galaxies at this redshift to determine how the galaxies in this 
object compare to other co-eval (proto-)clusters and to field 
galaxies. 

4.1. Color Evolution and Formation Epoch 

As discussed above, CIG J0218. 3-0510 is dominated by a 
strong red sequence of bright galaxies. These are intrinsically 
very luminous for their redshift, z = 1.62. The top axis of 
figure |6] compares the measured /-band magnitudes to the 
"characteristic" magnitude, J*, for the luminosity function of 
galaxies in the Coma cluster, evolve d passively backwards i n 
time to z = 1.62 using the model of Ide Propris et al.] (119991) . 
The brightest galaxies in CIG J0218. 3-0510 correspond to 
J* -I to /* - 1.5 mag. The descendants of these galaxies at 
a minimum will be super-L* cluster galaxies by z ~ 0, even 
without subsequent merging or star formation. 

The red galaxy colors imply that they contain older stel- 
lar populations. Figure |6] illustrates the expected color of 
a composite stellar population formed at z/ = 2.4 with a 
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Figure 6. Color-magnitude diagram of all galaxies within A < 2' of the cluster center, correspo nding to 1 Mpc at z = 1 .62. The top axis shows magnitudes in 
units of J*, the characteristic luminosity evolved from the measured values for the Coma cluster jde Propris et al.|[l999l) to ; = 1.62. The symbol size denotes 
that likelihood that the galaxy lies at the cluster redshift based on the integrated redshift probability, P,, as indicated in the legend. The smallest gray squares all 
have Vz < 0.3, and are likely foreground or background field galaxies, unassociated with the cluster Larger, red squares correspond to galaxies with V- > 0.3, 
and likely associated with the cluster The long-dashed line shows a fit to the red sequence (see § 4). The red dot-dashed line shows the expected color of a 
stellar population observed at z = 1.62 formed in a short burst at zf = 2.4 as described in the text. Circles denote those objects spectr oscop ically confirmed from 
1.62 < z< 1.65. These are bluer galaxies owing to the presence of [0 11] emission, facilitating their spectroscopic identification (see § 12.5) . The X symbol shows 
an object with a spectroscopic redshift outside this range. The short dashed lines indicate the 5<j magnitude limits derived for point sources in the imaging data. 
A strong red-sequence is apparent, which is dominated by galaxies with high Vz ■ 



Star-formation rate that evolves by an e-folding timescale 
r = 0.1 Gyr to z = 1-6 2 using the 2007 version of the 
iBruzual & Charloi (l2003h models. This model produces a 
(z-J) color similar to what we measure for the brighter red 
galaxies in this cluster (using the Bruzual & Chariot 2003 
models produces z-J colors within 0.02 mag for these ages 
and formation redshift). There is also no indication that 
CIG J0218. 3-0510 contains a population of blue, luminous 
galaxies. The brightest galaxies within a projected separation 
of r < 1 Mpc of CIG J0218.3-0510 with (z-J) < 1.5 mag and 
Vz > 0.3 have fainter magnitudes, typically sub-/*. Because 
the mass-to-light ratios of the stellar populations increase as 
they redden and fade over time, the dominant population of 
red-sequence galaxies could not form directly from the less- 
luminous blue galaxies. Interestingly, the number of red se- 
quence galaxies appears to decline toward the faint end of the 
red sequence, seemingly at magnitudes well above the detec- 
tion limit (see fi gure |6]i. This observation is similar to the 
findings of lTanaka et al.. (.2007.) . who observed a deficit of red 
galaxies around a z = 1 .24 cluster. These results imply a strong 
evolution in the faint end of the red-seque nce luminosity 
functio n, extending the relation measured by iRudnick et al.l 
(120091) for moderate redshift clusters to higher redshift. It also 
appears that there is a lack of faint blue galaxies with P. > 0.3 
but no such lack is seen in the V- < 0.3 galaxies, indicating 



that the photometry is indeed complete to the level indicated 
in the figure. Further studies of clusters at these redshifts and 
detailed modeling are needed to confirm this. 

The zeropoint, scatter, and slope of the red sequence 
itself provide constraints on the form ation timescales of 
the galaxies' stellar po pulations (e.g., iBower et al.l Il992t 
lAragon-Salamanca et al.|[T993i) . To compare the galaxies in 
CIG J0218. 3-0510 at z = 1.62 to those in lower redshift clus- 
ters, we converted the o bserved (z- J) co lors to rest-frame 
(U-B) colors following iMei et al.l (120091) . We find no evi- 
dence for evolution in the slope of the red sequence between 
CIG J0218. 3-0510 and lower redshift galaxy clus ters, con- 
sisten t with the results of other studies (see, e.g., IMei et all 
120091) . The top panel of figure|7]shows the rest-frame (U-B) 
color of CIG J0218.3-0510 and lower redshift clusters taken 
from the literature at a fixed rest-frame absolute magnitude, 
Mb = -21.4 mag. However, there is strong evolution in the 
rest-frame (U-B) colors, which become redder with decreas- 
ing redshift, consistent with passive evolution. The models in 
figure |7] illustrate the expected evolution for stellar popula- 
tions with different formation redshifts (zf = 2.0, 2.5, 3.0, and 
5.0), with a star-formation e-folding timescale, r = 0.1 Gyr, 
and solar metallicity. The (U-B) color at Mg = -21.4 mag 
for CIG J0218. 3-0510 is consistent with stellar populations 
formed between 2.0 < z < 2.5. 
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Figure 7. The top panel shows the evolution of the rest-frame U -B color 
measured at rest-frame Mg = -21.4 mag for red-sequence galaxies in clusters 
as a func tion of redsh ift. Smaller circles show the results from lower-redshift 
clusters IBower et al 1992; Ellis et al. 1997; van Dokkum et al. 1998a U 
1200(1 lMeretani2009l) . The large box point shows the value derived for the 
cluster CIG J0218. 3-0510 at 1.625 derived here. The curves show the ex- 
pected evolution of a stellar population with solar metallicity formed with an 
e-folding timescale t = 0. 1 Gyr at the formation redshift indicated. The bot- 
tom panel shows the evolution of the scatter in the rest-frame U—B color for 
red-sequence galaxies as a function of redshift. Symbols are the same as the 
top panel. The scatter in the colors of red-sequence galaxies increases with 
redshift. The right axis of the bottom panel indicates the expected scatter in 
the rest-frame (U -B) color for a red sequence observed at ; = 1.62 for a 
given mean luminosity-weighted formation redshift (see text). The scatter 
in the (U — B) color for the z = 162 cluster corresponds to a luminosity- 
weighted formation redshift oizf = 2.25 — 2.45, consistent with the evolution 
in the {U —B)o color. 



We model the scatter in the colors of the r ed-sequence 
alaxi es in CIG J0218. 3-0510 by following iHilton et al.l 
20091 and references therein). We construct a series of 
composite stellar populations with star-formation e-folding 
timescale of t = 0. 1 Gyr and solar metallicity. We tested other 
values for the metallicity (ranging from 0.2 Zq and 2.5 Zq), 
but found that they did not reproduce simultaneously the ze- 
ropoint and scatter in the rest-frame (U -B) colors for a con- 
sistent formation epoch. In our model each galaxy forms stars 
in a single burst of star formation that begins at an initial for- 
mation redshift, zf, and lasts for a duration At. We allow the 
duration Af to vary from Af = to a maximum equal to the 
lookback time from zf z= 1.62, the redshift at which the 
cluster is observed. We allow for a range of formation red- 
shift (with a maximum formation redshift taken to be when 
the lookback time was 4 Gyr), where at each, zf, we con- 
struct a simulated sample of 10^ galaxies with ages assigned 
randomly from a uniform distribution with < f < Af . From 
these simulated galaxies we compute the scatter in the sim- 
ulated color distributions and the luminosity-weighted age. 
Then, for a given measurement of the scatter in the colors, we 
infer the corresponding luminosity-weighted age and thus ob- 
tain an estimate for the formation redshift of the red-sequence 



galaxies in the cluster 

We calculated the scatter for the red sequence galaxies in 
CIG J0218.3-0510, including those galaxies with integrated 
redshift probability P- > 0.3 within a projected radius of 1 
Mpc (2 arcmin) of the cluster center We used an interative 
rejection algorithm to include only red galaxies within 2.5(7 
of the derived red sequence . We (ierive the sc atter about the 
median absolute deviation (iBeers et al.lll990l) . which yields 
a(U -B) = 0.136 ±0.024. We compare this to the red se- 
quences of other galaxy clusters and our models in the bot- 
tom panel of figure |2l There is a marked increase in the 
scatter in the {U -B) color with redshift, which is expected 
because the cluster galaxies are observed closer in time to 
when they formed their stellar populations (e.g., iMei et all 
120091: IHilton et all 120091) . The right axis of the panel in- 
dicates the modeled relationship between (t{U - B) and the 
formation redshift, zj for the model stellar populations ob- 
served at z = 1.62. The scatter in the rest-frame (U - B) 
color for CIG J0218. 3-0510 corresponds to range of forma- 
tion redshift, 2.25 £ z < 2.45, using the 2007 version of the 
iBruzual & Chariot (l2003h stellar population synthesis models 
(using the 2003 models has a negligible change on the forma- 
tion redshifts for the stellar population ages involved here). 
We note that this is an upper limit as we have made no at- 
tempt to remove the color measurement uncertainties from the 
intrinsic scatter. If we have overestimated the scatter, the for- 
mation redshift will be higher. Given the agreement between 
the formation epochs derived from the red-sequence zeropoint 
and scatter, we do not think this is a serious effect. 

Therefore, both the zeropoint and scatter of the rest-frame 
(U - B) colors for red-sequence galaxies in CIG J0218.3- 
0510 imply a formation epoch of z/ — 2.25-2.45 (a look- 
back time of 1.0-1.3 Gyr from z = 1.62). This corresponds 
to the last major star-formation episode in these red cluster 



galaxies. This is a similar formation epoch as found in many 

;r galaxies (e.g.. iMei et alj|2009l : 

IHilton et al J 120091: iRosatTetallllOOgl) . 



studies of other z > 1 cluster galaxie s (e.g.. lMei et al 

although at least one 



cluster (XMM J2235. 3-2557 at z = 1.39) has g alaxies with 
formation epochs ranging from z/'--"2toz/~6 ( iRosati et alj 
I2009h . The formation epoch of CIG J0218. 3-0510 is also 
consistent with the colors and spectra l indices of galaxies 
in lower-redshift massive clusters (e .g., IStanford et al]ll998l : 
Ivan Dokkum & van der Mareil 120071) . Furthermore, studies 
show that a high fraction of massive galaxies at z > 2 
have high levels of star formation (e.g.. iKriek et al.l 120061 : 
iPapovich et al.ll2006h . and these are likely consistent with the 
evolution of the cluster galaxies. These properties all support 
the conclusion that the galaxies of CIG J0218. 3-0510 will be- 
come typical galaxies found in rich clusters at later times. 



4.2. Dynamical Mass Estimate 

As discussed above, it is unclear if the cluster galaxies of 
CIG J0218. 3-0510 sample a virialized, relaxed dark-matter 
halo, or if this object is in the process of assembling through 
mergers. Under the assumption that the galaxies are virial- 
ized, we use the velocity dispersion to provide a crude esti- 
mate for the dynamical mass of CIG J0218. 3-0510 a s a r efer- 
ence, although we qualify this with the caveats in § 12.51 For 
a velocity dispersion, a = 860 ± 490 km s"' (see § 12.5b . the 
corresponding dynamical mass estimate is Mdyn ^ 3rvii-o'" ~ 
4 X 10'"* Mq for a virial radius of 0.9 Mpc. We derive the same 
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dynamical mass estimate using M200 (ICarlberg et al.|[T997l)n 
Even so, we caution the reader that the dynamical mass may 
be highly uncertain given the statistical errors and systematic 
biases. 

This mass is a factor of ^4 large r than th e limiting halo 
mass of the IRAC-selected clusters (iPapovich 2008). and if 
accurate then it implies that CIG J0218. 3-0510 resides in one 
of the most overdense environments in the SWIRE survey. In 
this case, based on th e expected growth of dark-matter haloes 
(ISpringeletal.ll200l the mass of CIG J0218.3-0510 should 
increase to at least lO'*" Mq by z = 0.2, becoming a rich cluster 
of galaxies comparable to the Coma cluster. 

The high dynamical mass estimate of CIG J0218. 3-0510 
corresponds to a large gravitational potential well. If this is 
the case, then the hot diffuse gas of the ICM should emit 
significant X-ray luminosity. This is supported by the X- 
ray detection we infer at the location of CIG J0218. 3-0510 
in § 12.61 Assuming that all the detected X-ray emission 
originates from the ICM we infer a mean temperature of 
~3 keV and virial mass of 1.1 X 10"* Mq based on the lo- 
cal luminosity-tempera t ure and mass-temp erature relations 
jArnaud & Evrardlll999l:P^rnaud et al.l l2005!). This is roughly 
a factor of 4 lower than the dynamical mass derived above, 
but within the large error budget. While the X-ray detec- 
tion of CIG J0218. 3-0510 supports the high derived virial 
mass, obtaining much deeper X-ray data is required for a self- 
contained analysis of the ICM properties and mass determina- 
tion. 

5. SUMMARY 

We report the discovery of a spectroscopically confirmed 
galaxy cluster at z=1.62 located in the SWIRE XMM-LSS 
field. This cluster candidate was selected solely as an over- 
density of sources with red Spitzer/IRAC colors, satisfying 
([3.6]- [4.5])ab > -0.1 mag, with no other selection criteria 
imposed. Photometric redshifts derived from SXDS {BViz- 
bands), UKIDSS-UDS (//iT-bands), and SpUDS (3.6-8.0 mi- 
cron) for the galaxies in and around this cluster show that this 
structure corresponds to a galaxy surface density of sources 
at z = 1.6 that is > IQa times the mean surface density at this 
redshift. Furthermore, our recent analysis of existing XMM 
data on this cluster provides a weak but unambiguous detec- 
tion compatible with the expected thermal emission from such 
a cluster This confirms that our selection of overdensities 
of sources with red IRAC colors identifies galaxy clusters at 
z> 1.3. 

We obtained spectroscopic observations of galaxies in the 
cluster region using IMACS on the Magellan telescope. We 
measured redshifts for five galaxies in the range z= 1.62-1. 63, 
all within 1.4 arcmin (< 0.7 Mpc) of the cluster center. In ad- 
dition, we measured spectroscopic redshifts for two sources 
with z=1.62-1.63 within 1.4-2.8 arcmin (0.7-1.4 Mpc). 
The cluster appears to be dominated by red galaxies, with 
(z-J) > 1.7 mag. The photometric redshift distributions for 
the brightest red galaxies are centrally peaked at z = 1 .62, co- 
incident with the spectroscopically confirmed galaxies. 

The z-J versus J color-magnitude diagram of the galaxies 
in this cluster shows a strong red-sequence, which includes 
the dominant population of red galaxies. The rest-frame 
(U - B) color and scatter of galaxies on the red-sequence are 

Here Miqo is the mass of the dark matter halo within a spheiical volume 
defined by radius rioo where the average density is 200 times the critical 
density at the observed redshift. 



consistent with a mean luminosity-weighted age of 1.0- 1.3 
Gyr, yielding a formation redshift z/ = 2.35 ± 0. 10, and corre- 
sponding to the last major episode of star formation in these 
galaxies. 

We provide a crude estimate of the dynamical mass of 
CIG J02 18.3-05 10, although this result is highly uncertain due 
to systematics in the data and owing to the unknown dynam- 
ical state of CIG J0218. 3-0510. Under the assumption that 
the dark matter halo of CIG J02 18. 3-05 10 is virialized and re- 
laxed, we estimate a dynamical mass based on the measured 
velocity dispersion, M^ya ^ 4 x 10'"* Mq within a virial radius 
of rvii ^ 0.9 Mpc. Our reanalysis of the XMM X-ray emission 
of this cluster favors a lower virial mass, 1.1 x 10'"* Mq, but 
consistent within the uncertaintites. If these masses are ac- 
curate, then we expect CIG J0218. 3-0510 to evolve to a rich 
cluster with M ^ 10'^ Mq at z = 0.2 similar to the Coma clus- 
ter Further spectroscopic and multiwavelength observations 
of galaxies in CIG J0218. 3-0510 are needed to constrain bet- 
ter the mass measurement. 
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